Transition from spatial coherence to incoherence in coupled chaotic systems 
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We investigate the spatio-temporal dynamics of coupled chaotic systems with nonlocal interac- 
tions, where each element is coupled to its nearest neighbors within a finite range. Depending upon 
the coupling strength and coupling radius, we find characteristic spatial patterns such as wave-like 
profiles and study the transition from coherence to incoherence leading to spatial chaos. We analyze 
the origin of this transition based on numerical simulations and support the results by theoretical 
derivations identifying a critical coupling strength and a scaling relation of the coherent profiles. 
To demonstrate the universality of our findings we consider time-discrete as well as time-continuous 
chaotic models realized as logistic map and Rossler or Lorenz system, respectively. Thereby we 
establish the coherence-incoherence transition in networks of coupled identical oscillators. 

PACS numbers: 05.45.Xt, 05.45.Ra, 89.75.-k 

Keywords: nonlinear systems, dynamical networks, coherence, spatial chaos 



I. INTRODUCTION 

While isolated nonlinear systems are fairly well under- 
stood, much less is known about emergent dynamics of 
coupled systems [THS]. Complete synchronization, also 
known as zero-lag synchronization, where all nodes of 
a network behave identically, is just the simplest sce- 
nario. Nevertheless, it has important applications, for 
instance, as chaos synchronization for secure communi- 
cations [3 IH] . Other types of synchronization have been 
found as well, including phase-lag, cluster, or generalized 
synchronization jj. Recently a new form of dynamics 
called chimera state triggered a lot of research [5l419j. 
These chimera states arise in regular networks of identical 
phase oscillators and consist of parts with high spatial co- 
herence and regions where spatial coherence is lost. They 
can be realized with carefully prepared initial conditions 
and a nonlocal network topology, which both arc crucial 
ingredients for these states. Similar phenomena, local- 
ized patterns surrounded by a chaotic background, were 
found in populations of threshold elements with nonlocal 
interaction, related to a continuity-discontinuity transi- 
tion pni- 

In this paper we investigate in detail the transition 
from spatial coherence to incoherence in both time- 
discrete and time-continuous chaotic systems, using the 
well-known models of the logistic map and the Rossler 
system, respectively. See also Ref. [H] for a first report 
on this transition. We consider a ring topology and use 
the coupling strength and coupling range as bifurcation 
parameters. At the transition we observe partially in- 
coherent chimera-like states that we identify as spatial 
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chaos da US]. 

This paper is structured as follows: Sect ion [II] intro- 
duces the model equations. Sections |III| and |IV| describe 
the general coherence-incoherence scenario and numer- 
ical results on the coherence-incoherence transition, re- 
spectively. Analytical results are presented in Sec. |V] In 
Sec. IVII we demonstrate that the observed scenario can 
also be found in the Lorenz model. Finally we conclude 
in Sec. |VlIl 



II. MODELS 

In order to investigate the universality of the 
coherence-incoherence transition for N coupled oscilla- 
tors we consider different nonlinear models: As chaotic 
systems we employ the time-discrete, logistic map as well 
as the Rossler system. The latter is paradigmatic for 
chaos described by ordinary differential equations. 

As a control parameter we use the coupling strength 
a and the coupling radius r — P/N, where P denotes 
the number of nearest neighbors in each direction on a 
ring. Varying the value of P from 1 to N/2 we cover the 
range of regular network topologies from local (nearest- 
neighbor) to global (all-to-all) coupling. 



A. Logistic map 

As an example for time-discrete systems, we study cou- 
pled maps: 

4^' = fi4) + ^ E [/ (4) (1) 

J=i-P 



2 



where Zi (i = 1,...,7V) are real dynamic variables and 
t denotes the discrete time. In the coupling term the 
index i is periodic mod N. The function /(z) is a local 
one-dimensional map that we choose as the logistic map 
f{z) = az{l — z). We fix the bifurcation parameter a at 
the value a — 3.8 unless stated otherwise. This is a value 
that yields chaotic behavior with a positive Lyapunov 
exponent A « 0.431. 



B. Rossler system 

In order to show that the coherence-incoherence tran- 
sition scenario is universal, we will also discuss nonlo- 
cally coupled networks of time-continuous systems. As 
a prominent example for chaotic systems, we consider 
the Rossler model. In analogy to Eq. ([!]) the compound 
system is given by 
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FIG. 1. (Color online) Blowout bifurcation curves in the 
(r, (T)-parameter plane for A'^ = 100 coupled logistic maps. 
The system parameter is chosen as a = 4.0, 3.8, and 3.6 for 
the solid, dashed (red), and dotted (blue) curve, respectively. 
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Zi = b + Zi{xi - c) -I- — i^J - ^i) ; 

j=i-P 



with i = 1,...,A''. Again the index i is periodic mod 
N to realize the underlying ring topology. The system's 
parameters are chosen as a = 0.42, b — 2, and c = 4 
showing chaotic behavior in the uncoupled case. 

The 3A^ coupled Eqs. ^ can also be written in com- 
pact form as follows 



F(X) 



2P 



(G®H)X 



(3) 



with X — {xi,yi, zi, . . . ,XN,yN, zn) and a 3 x 3 cou- 
pling matrix H that determines, which components are 
transmitted. The N x N matrix G describes the network 
topology. In the case of Eqs. we have the identity 
matrix for H, and a circulant matrix with cyclically per- 
mutated rows (— 2P, 1, . . . , 1, 0, . . . , 0, 1, . . . , 1) is chosen 

p p 
for G such that gu — ~2P. 

This notation has been employed in the framework of 
the master stability function in order to investigate com- 
plete, potentially chaotic synchronisation Based on 
a variational equations this formalism allows for calcula- 
tion of the Lyapunov exponent transverse to the synchro- 
nisation manifold |251 126j . The master stability function 
technique has been used in various contexts ranging from 
generic models [171 HH] to coupled lasers [IHl [SD] and neu- 
roscience [3Tti33] . It can also be extended to cluster syn- 
chronization |34| , but does not give conclusive results for 
the purpose of the investigation of spatially incoherent 
states. 



III. STRUCTURE OF COHERENCE TONGUES 

Let us discuss first the coherence-incoherence transi- 
tion for coupled logistic maps described by Eq. ([T]). For 
large values of coupling radius r and coupling strength a 
we observe complete chaotic synchronization. There the 
oscillators behave identically {zsync = zi = ... = zj^), 
but chaotically in time following the dynamics of /. 
Using the fact that the coupling matrix of system ([T]) 
is a circulant matrix, with cyclically permutated rows 
(— 2P, 1, 1, 0, 0, 1, 1), we obtain analytically the 

p p 
values for — 1 transversal Lyapunov exponents of the 
completely synchronized solution Zgy^c'. 
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m = 1, N — 1, where Ay > denotes the longi- 
tudinal Lyapunov exponent corresponding to the one- 
dimensional local map, and the second term on the right- 
hand side is derived from the eigenvalues of the Jacobian 
matrix of the system, using the fact that coupling matrix 
is circulant. 

With a decrease of the coupling radius r or the coupling 
strength a the completely synchronized solution loses its 
stability in a blowout bifurcation, where the largest trans- 
verse Lyapunov exponent becomes positive. The position 
of the blowout bifurcation line in the (r, (T)-plane depends 
on the local dynamics of the system elements. Figure [T] 
shows examples for the blowout bifurcation curves for 
a — 3.6, 3.8, and 4.0. For any of these values the local 
logistic maps operate in the chaotic regime. The region 
of complete chaotic synchronization is the region above 
the blowout bifurcation curves. 
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The completely synchronized state represents the sim- 
plest example of spatial coherence. To be precise, we will 
call a network state z*, i — I, ■■■,N, coherent on the ring 
as TV — >■ oo, if for any point a; G 5^ the following limit 
holds: 



lim lim 

TV-i-oo t-i-oo 



sup 



0, for (5^0, (4) 



where {x) = {j : < j < N, \j/N - x\ < 6} denotes 
a network neighborhood of the point x. Intuitively it 
means that the profile becomes smooth in the thermo- 
dynamic limit N ^ oo. If the limit (4| does not vanish 
for i5 — 0, at least for one point x, the network state is 
considered incoherent. 

Below the blowout bifurcation curve, several further 
regions of coherence appear as the control parameters a 
or r decrease. Figure [2] depicts these regions as shaded 
(color) tongues for the logistic map and the Rossler sys- 
tem in panels (a) and (b), respectively. The blowout 
bifurcation line is added as dashed curve (BB), and the 
regime of complete, chaotic synchronisation is indicated 
by the light hatching. 

A coherent state has a smooth profile characterized by 
the number of maxima, i.e., the wave number k. Only 
regions for wave numbers k = 1, 2, and 3 are shown. The 
corresponding coherent states zj are shown as snapshots 
in the insets of Fig.jSJa) for the logistic map. Further de- 
crease of r yields additional thin higher-order regions fol- 
lowing a period-adding cascade fc = 4, 5, . . . (not shown). 
Inside the regions, the states are coherent in space and 
periodic in time and undergo a period-doubling cascade 
in time as r or cr decrease. In the parameter space be- 
tween the coherence regions the network dynamics re- 
mains coherent but not periodic anymore. The states 
alternate chaotically between the adjacent fe-states and 
thus exhibit chaotic itineracy [2j |35] . The combination 
of period-adding in space and period-doubling in time 
represents a remarkable feature of networks of coupled 
chaotic oscillators with nonlocal coupling. 

Figure shows stability diagram for the coupled 
Rossler systems again in the (r, cr)-plane. It is easy to 
see that the overall structure of the diagram is similar 
to the one obtained for nonlocally coupled logistic maps 
(Fig.[2ja)). In analogy to the maps, there exist a blowout 
bifurcation boundary (BB curve) and a sequence of co- 
herence tongues with wave numbers k indicating the uni- 
versality of the general dynamical scenario. In the shaded 
(color) regions shown in Fig.|2jb) the coherent states with 
regular (periodic or quasiperiodic) dynamics are the only 
stable solution. The bright (yellow) regions inside the co- 
herence tongues below the stable period-2 tongues refer 
to multistability and more complex behavior including 
coexisting period-4 and torus solutions in between the 
period-2 and period-4 regions. The dark (blue) region 
is the region of stability for the solutions with periods 8 
and larger. Since these regions appear to be very thin, 
we combine them for better visualization. 

For further insight of the transitions within the co- 
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FIG. 2. (Color online) Coherence regions in the (r, a) param- 
eter plane for A'^ = 100 logistic maps (a) and Rossler systems 
(b), labeled by the wave number k. Gray scale (color code) in- 
side the coherence regions distinguishes different periods and 
the coherence-incoherence bifurcation (CIB) curve separates 
regimes with coherent and incoherent dynamics. The light 
hatched region bounded by the blowout bifurcation curve BB 
refers to completely synchronized chaotic states. The insets 
in panel (a) display snapshots of typical coherent states. Sys- 
tem parameters: panel (a) a = 3.8; panel (b) a — 0.42, b — 2, 
and c = 4.0. The vertical dashed lines refer to values of a 
in Fig. [3] Triangles denote parameter values used to describe 
the bifurcation scenario in Figures [4] and |5] 



herence tongues Figs. |3]^a) and (b) depict cuts through 
the (r, (T)-parameter space for a fixed coupling radius 
r — 0.12 and 0.24, respectively (see dashed lines in 
Fig. [2jb)). The black and gray (red) values refer to the 
variable j/i at a Poincare section (xi = 0) for a very slow 
decrease and increase of the coupling strength a with 
steps A(T = 0.001, respectively. For each value of cr 300 
iterations are shown, where each initial condition is cho- 
sen as the final configuration of the previous value of cr. 
Panel (a) shows a sequence of chaotic, period-4, torus, 
period-2, coexisting period-2 and period-4, period-4, and 
again chaotic dynamics. In panel (b) there appears an 
additional regime of chaos in between the period-2 and 
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FIG. 3. (Color online) Bifurcation scenario at cuts through 
the fe = 1 coherence tongue of Fig. [2|b) at (a) r = 0.12 and 
(b) r = 0.24. The Poincare section {xi = 0) of yi is shown for 
the Rossler systems in dependence on the coupling strength 
a. The black and gray (red) points correspond to a slow down 
and up sweep of a, respectively. Long transients up to 3000 
oscillation periods have been discarded. Other parameters as 
in Fig. [2] 



period-4 region. 

The coupling radius appears to be the crucial parame- 
ter in the system with nonlocal couplings. With increas- 
ing the system size N and a proportional increase of the 
number of coupled neighbors, we obtain similar coherent 
solutions. If with increasing TV the number of coupled 
neighbors P does not change, then the coupling radius 
decreases, and the point of operation (r, a) moves in the 
stability diagram (Fig. [2]) to the left, obtaining coherent 
solutions with larger wavme numbers or itineracy regions 
between them. 

The bifurcation diagrams presented in Fig.[2]have been 
obtained by means of numerical simulations. Note, how- 
ever, that in the thermodynamic limit A'^ — ^ oo it is also 
possible to obtain some analytic results. Details on such 
derivations are given in Sec.|Vj e.g., by identifying a scal- 
ing relation of the coherent profiles. 



IV. COHERENCE-INCOHERENCE 
TRANSITION 



The described combination of period-adding in space 
and period-doubling in time for coherent solutions is 
valid for sufficiently large coupling strength (cr > 0.4). 
In Ref. j21J we have shown that with decreasing cou- 
pling strength a a coherence-incoherence bifurcation is 
observed in the time-discrete system ([T]). The typical sce- 
nario of this transition is shown in Fig. [4] for the logistic 
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FIG. 4. Snapshots of the dynamics of the logistic map for 
different coupling strengths a = 0.43,0.4,0.32,0.3,0.2, and 
0.1. The gray regions labeled ai and a2 refer to multistable 
incoherent parts of the profile. The coupling radius is fixed 
at r = 0.32. Other parameters as in Fig. [2] 



map: The smooth wave-like solution profile z* splits up 
into upper and lower branches. Furthermore, two nar- 
row boundary layers of incoherence appear around the 
points xi and X2 in Fig. |4|^b). They are indicated by 
shaded stripes ai and a2 in Figs.gc) and (d). The width 
of the incoherence stripes increases for smaller coupling 
strength a (see Fig. |4]^d)) until the dynamics becomes 
completely incoherent as depicted in Figs. |4|^e) and (f). 

Below the bifurcation parameter value a « 0.40 nu- 
merous hybrid, that is, partially coherent states arise. 
They are characterized by coherence on some intervals of 
the ring and incoherence on the complementary in- 
tervals. Figs.|4][c) and (d) display exemplary realizations 
of these partially coherent states. Within the incoher- 
ence intervals, i.e., ai and a2, any combinations of the 
upper and lower states are admissible for appropriately 
prepared initial conditions. 

Thus the coherence-incoherence bifurcation leads to 
spatial chaos ,22. .23j that appears first at narrow inco- 
herence intervals and, with decreasing coupling strength, 
encompasses the whole ring. Therefore we identify a 
chimcra-like state of coexisting coherent and incoher- 
ent regions as a transitional state in the coherence- 
incoherence bifurcation scenario. Opposed to previ- 
ously reported chimera states in time-continuous systems 
[21 [ini however, the temporal behavior in the 

present case is periodic rather than chaotic. The com- 
plexity arises due to the huge variety of multistable inco- 
herent states and corresponds to permutations of the se- 
quence of upper and lower local states. With further de- 
crease of the coupling strength, the chimera states cease 
to exist and we find completely incoherent behavior. 

The findings above are based on numerical simulations 
of coupled chaotic logistic maps. Even if the local logistic 
maps operate in the periodic regime, one can observe a 
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sequence similar to Fig. |4] as reported in Ref . [21] • Note 
that, when the local dynamics of the system is periodic, 
the region of stability of the completely synchronized so- 
lution covers almost all of the (r, cr) parameter plane, and 
the system is characterized by multistability. 




FIG. 5. (Color online) Snapshots of the variables Xi (upper 
panels, X50 as black dot) and space-time plots (lower panels) 
for the Rossler system. The coupling strength is chosen as 
a = 0.12,0.1,0.09, and 0.06 in (a) to (d), respectively. The 
coupling radius is fixed at r = 0.28. Other parameters as in 
Fig. [2] 

Figures [5] and [6] display the typical scenario of 
coherence-incoherence transition for nonlocally coupled 
Rossler systems, similar to the one described above for 
the logistic map. In the upper panels of Fig.jsj snapshots 
of the variables shown. Here all nodes are shown 

in gray, and an exemplary node (i = 50) is highlighted as 
a black point. Similarly to the coupled maps, the smooth 
profile of the coherent (fc = 1) solution of Fig.jsjja) breaks 
up into two clusters in Fig. [5][b). With further decrease 
of the coupling strength a partially coherent solution 
appears in the system as shown in Fig. [sjc). Finally 
we observe a chimera-like state around the upper clus- 
ter state in Fig. [Sjjd) , which combines spatially coherent 
and spatially chaotic parts. The lower panels of Fig. [5] 
show space-time plots. Here the spatial structure and the 
breaking up of the profile is clearly visible. Furthermore 
one can already infer the temporal behavior that will be 
discussed in the following. 

Figure |6] shows three-dimensional phase portraits of 
the Xi variables and the superimposed time series of all 



Xi{t), i = 1,...,N of the Rossler systems for the same 
parameter values as in Fig. [5] Again the time series of 
node i = 50 is highlighted as black (red). Note that 
in Figs, and (c) the time series of the points from 
the incoherent regions are grouped into two clusters that 
are formed by the two coherent parts of the solution. In 
Fig. [6]Jd) one can visualize the chaotic dynamics already 
in the short time interval shown. In the time domain one 
can observe a period-doubling bifurcation in Figs. [6]^ a) to 
^c) eventually leading to chaos in Fig. ^d) (lower pan- 
els). At the same time the spatial irregularity becomes 
more pronounced as described above. 



V. ANALYTICAL RESULTS 

After the numerical exploration of the coherence- 
incoherence transition, we provide some analytical 
derivations in this Section. We will discuss an approxi- 
mation for the value of the critical coupling strength at 
the coherence-incoherence transition and present a scal- 
ing of the profiles in different coherence tongues, which 
holds in the limit A^ — > 00. 



A. Critical coupling strength 

In the following, we derive an approximation of the 
critical coupling strength at the transition from coherent 
to incoherent states, see the line CIB in Fig.[2];a). We will 
restrict the investigation to the logistic map introduced 



in Sec. II A but a similar argument holds for the Rossler 
system as well. 

From a geometrical point of view and in the thermo- 
dynamic limit N ^ 00, coherent solutions z* approach a 
smooth profile z^{x) of the spatially continuous version 
of Eq. given by 

fX-\-r 



z'+'{x)^f{z\x)) + 



2r 



[f {z\y))~f{z\x))]dy. 

(5) 

A transition from coherence to incoherence occurs when 
the respective solution profile z*{x) becomes discontinu- 
ous in some points x of the ring S^. 

Let us now discuss in more detail this coherence- 
incoherence bifurcation. Consider a solution of sys- 
tem ([5| with wave number k = 1 and period-2 dynamics 
in time. Hence we can reduce the dynamics by even and 
odd time steps zo{x) and zi{x), respectively. This leads 
to 



x-\-r 

_^(x)^{l-a)f{z,{x)) + ^ J f{z,{y))dy 



(6) 



with j = 0, 1. Taking the spatial derivative yields 

z[_j{x)={l-a)f {z,{x))zr{x) 

+ Y^[f{zA^ + r))~.f{z,{x~r))]. (7) 
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FIG. 6. (Color online) Phase space trajectories (upper panels) and time series Xi{t), i = 1,2, ...,A'^, (lower panels) for the 
Rossler system. The coupling strength is chosen as cr = 0.12, 0.1, 0.09, and 0.06 in panels (a) to (d), respectively, corresponding 
to Fig. [5] The exemplary time series xso{t) is highlighted in black (red). The coupling radius is fixed at r = 0.28. Other 
parameters as in Fig. [2] 



At the point x where the smooth profile breaks up, 
the spatial derivative becomes infinite. Considering that 
Zq{x), z[{x) diverge to infinity, we can neglect the cou- 
pling term on the right-hand side of Eq. ([t]). The main 
contribution comes from the first term. 

Multiplying the equations for even and odd time steps 
we obtain: 

z',{x)z[{x) = [(1 - aff (zoix)) f z',{x)z[{x), 

(8) 

which yields the following condition 

l = {l-aff{zo{x))r{z,{x)). (9) 

Keeping in mind that the local dynamics is governed by 
the logistic map f{z) — az(l — z) and that its derivative 
is equal to f (z) = a(l — 2z), we introduce the function 

G{x) = (1 - afa^ (1 - 1zq{x)) (1 - 2zi(a;)) - 1 (10) 



as a deviation from condition (|9]). 

For fixed system parameters r and a that correspond 
to the region of coherent solutions with periodic behav- 
ior of period 2 in time, and thus fixing zq(x) and z\(x\ 
we can calculate G(x) and determine values x, where 
G[x) = 0. The obtained values x are the locations at 
which the smooth profile break up into two parts. Fig- 
ure [t] illustrates G{x) for a fixed coupling radius r = 0.32 
for the finite case of 100 elements. Numerical evaluation 
shows that for a coupling strength a « 0.41 the function 
G approaches zero in two points. Figure[7]^b) displays the 
corresponding snapshot of the solution. One can see that 
the profile breaks up at positions at which G vanishes as 
indicated by the vertical dashed lines. 

To derive an approximation for this coupling strength 
let z* be the fixed point of the local logistic map: z* — 
f{z*) = az*{l — z*), hence z* = 1 — 1/a. Here, for 
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FIG. 7. (a) Deviation from the incoherence condition Q for 
coupled logistic maps for a — 0.42 (dashed), 0.40919 (solid 
curve), G — 0.40 (dotted), (b) Snapshots for even (open dots) 
and odd (filled dots) solutions for cr = 0.40919 and r = 0.32. 
Other parameters as in Fig. [2] 



a — 3.8 we have z* « 0.737. Under the assumption 
zq{x) — z\(x) = z* if G = 0, we obtain an approximation 
for a: 



G(x) = [a(l-(T)(l-2z*)]' 



1 = 



(11) 
(12) 



convenience, we consider the case of time-continuous sys- 
tems, but our findings are also valid for maps as we will 
show later (Fig. |8|. 

The network model of a ring of N linearly coupled 
oscillators, each coupled to its P nearest neighbors. 



f(x,) 



cr 

2P 



^ (Hx,-HxO 



(13) 



becomes, in the limit iV — >■ oo, 

fU-\-T 



x(m) 



f (x(m)) + ^ \ [Hx(w) - Hx(u)] dw (14) 

^"^ Ju~r 



with coupling range r and coupling strength a as used 
in the previous sections, u e [0, 1] denotes the spatial 
position on the ring. Compare Eq. ([s]) for the spatially 
continuous version of time-discrete maps. The coupling 
scheme is determined by the matrix H. In the case of the 
Rossler system introduced by Eqs. ([2]), H is the identity 
matrix and the vector describing each node is given by 
Xi = {xi,yi,Zi). Considering the logistic maps of Eq. ([I]), 
the coupling term Hxj changes to f{zj), but the follow- 
ing argument still holds (Fig. [8]ja),(b)). 

In short, we will show that for fixed a the spatially 
fci-periodic dynamics of an oscillator at the position u 
with coupling range ri is identical to that at a rescaled 
position uki / k2 of the corresponding fc2-periodic state in 
a network with coupling range r2 = riki/k2. From this 
we can conjecture to find a fc2-periodic state at r2, if a 
stable fci-periodic state is found at ri, such that fc^rj = 
const. 



For a = 3.8 we get a « 0.44. 

From the numerical experiments we obtain that the 
transition from coherence to incoherence in system ([T]) 
occurs for a coupling strength cr close to 0.4, which 
roughly agrees with our approximation. The deviation is 
due to the finite number of nodes in the simulations. The 
derivation above was presented for wavenumber k = 1. 
Note that coherent solutions with higher wave numbers 
undergo the same bifurcation scenario, with the smooth 
profile breaking up into 2k parts. 



B. Scaling of coherent profiles 

As a second analytic result we will derive in the follow- 
ing a scaling relation that allows a mapping of coherent 
profiles in different regimes of coherence. These regimes 
were depicted in Fig. [2] by shaded (color coded) tongues 
labeled fc = 1, 2, . . . . Our considerations apply to rings of 
oscillators in the thermodynamic limit N ^ oo, N being 
the number of coupled oscillators. A comparison with nu- 
merical simulation shows that for sufficiently large finite 
iV the result holds in good approximation. For notational 
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FIG. 8. (Color online) Scaling of coherent profiles: Black 
points are snapshots of numerical simulations (a),(b) for cou- 
pled logistic maps with A'^ — 100, a = 0.7, r — 0.22 and 
r — 0.11, respectively, and (c),(d) for Rossler systems with 
N = 100, a = 0.6, r = 0.12 and r = 0.06, respectively. Gray 
(red) points in panels (b),(d) show rescaled k = 1 solutions 
from panels (a) and (c). Other parameters as in Fig. [2] 



Assume that a state Xs(m) with spatial period ki exists. 
For simplicity and notational convenience consider ki = 



8 



1. This state evolves according to 
±s{u) = HM^)) + ^ I ^ [Hx,(w)-Hx,(u)]dt;. 

Ju—r 

(15) 

In order to compare another fc-periodic state xj, with the 
state Xs, in the way described above, we assume 



yis{u) = Xfc 



u + m 



(16) 



with m = 0, 1, 2, . . . , fc — 1. Vice versa, scaUng the spatial 
coordinate of Xfe(u) with u G [0,1], as well as r, by a 
factor leads to 



±kik-'u) (xfc(fc-iu)) 

ka rk-\u+r) 



fc-i(«-r) 



[Hxfc(«)-Hxfe(fc-iM)] dv. 

(17) 

Substituting v = k^^w gives 
±k{k^^u) =f (xfc(fc~-^u)) 

(18) 



and, using Eq. (16), results in 



Xfc(fc ^U) = f (Xs(u)) + 



(T 

2^ 



u-\-r 



[Hxs(u;) — Hxs(ii)] d-Ti;. 

(19) 



Comparing the right-hand side of Eq (19) to Eq. (15 1, 
we find Xfe(fc~^w) = Xs(u), if the coupling radius r is 
appropriately rescaled by k~^. 

This derivation shows that one can recover the cou- 
pling parameters leading to a state with wave number ^2 
from the parameters of a ki state. This finding is illus- 
trated in Fig. [H] for ki — 1 and ^2 = 2 by corresponding 
snapshots of 100 logistic maps and 100 Rossler system 
in panels (a),(b) and (c),(d), respectively. Already for 
N = 100 there is strong agreement between the numeri- 
cally obtained and scaled profiles. In the case of logistic 
maps (Rossler systems), we choose the coupling strength 
as cr = 0.7 (0.6) and the coupling radius as r = 0.22 
(0.12) in Fig. jsf^a) (Fig. [8](c)) and, according to the scal- 
ing relation, as r — 0.11 (0.06) for the k2 — 2 profile in 
Fig. Mh) (Fig. ^d)). The gray (red) dots in Figs. Mb) 
and (d) depict the scaled profile of Fig. [8](a) and (b), 
respectively. 



VI. LORENZ MODEL 

To investigate how universal the results obtained so 
far are, we provide one more example of a chaotic 



(a) 8fo=8 




FIG. 9. Snapshots of nonlocally coupled Lorenz systems 1 20 1 
for different values of coupling strength a = 8, 6, 5.5, and 1 in 
panels (a) to (d), respectively. The coupling radius is fixed at 
r = 0.3 and the number of nodes is = 100. Parameters of 
the Lorenz system: ? = 10, p = 28, /3 — 8/3. 



time-continuous system that exhibits the coherence- 
incoherence bifurcation scenario described above. Fig- 
ure [9] displays the coherence-incoherence transition for 
nonlocally coupled Lorenz systems: 



^(y» - 


Xi) 4- 


a 

2P 


i+P 

j=i-p 


Xi) , 










iVj - 


Zi) - 


Vt 4 


-— E 




4 


a 

' 2P 


i+P 
J=i-P 


- Zi) , 



(20) 



where the local parameters are fixed at the values <j = 10, 
p = 28, (3 — 8/3. Fig. [9] reveals a striking similar- 
ity with the sequences of snapshots for nonlocally cou- 
pled chaotic logistic maps (Fig. |4| and Rossler systems 
(Fig. [5]). For fixed coupling radius and with decreas- 
ing coupling strength the smooth profile breaks up into 
two parts (Fig. |9][b)). Further decrease of the coupling 
strength leads to the appearance and growth of inco- 
herent parts, repeating again the typical scenario of the 
coherence-incoherence bifurcation described above. 



Additional numerical evidence shows that system ( 20 ) 
is characterized by high multistability involving different 
types of solutions. There are no clearly separated tongues 
in the (r, cr)-parameter plane present as in logistic maps 



9 



and Rossler systems, but there is strong overlap of the 
regimes of different wave numbers A; = 1, 2, 3, . . . and the 
region of spatially coherent synchronized chaotic dynam- 
ics (fc = 0). This is due to the completely different local 
dynamics, which is characterized by a double scroll and 
two symmetric saddle-foci. The coupling induces stable 
spatially coherent steady states in the form of wavelike 
profiles with wave numbers = 1, 2, 3, . . . (Fig. ^a.)). 

VII. CONCLUSION 

By both numerical and analytical means we have in- 
vestigated the transition from coherence to incoherence 
in coupled chaotic systems. In dependence upon the 
coupling strength and coupling range we have observed 
a combination of period-adding in space and period- 
doubling in time for nonlocally coupled chaotic maps and 
Rossler systems. An inspection of the spatial profiles 
reveals a coherence-incoherence bifurcation, where the 
smooth coherent profile of the solution breaks up into 
parts and gives rise to partially coherent solutions and 
chimera states. Upon further decrease of the coupling 



strength completely incoherent states are found as a sign 
of fully developed spatial chaos. A similar transition from 
coherence to incoherence is found for nonlocally coupled 
Lorenz systems. Our findings hold for both time-discrete 
and time-continuous systems indicating the universality 
of the transition. 

In conclusion, we have analyzed in detail the transition 
between spatial coherence and incoherence that appears 
in different systems with nonlocal couplings, and we find 
similar bifurcation scenario for nonlocally coupled logistic 
maps, Rossler and Lorenz systems. Beyond numerical in- 
vestigations, wc have determined analytically the critical 
coupling strength for this bifurcation, and established a 
scaling relation for profiles with different wave numbers. 
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